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A combined analyses data on isoscalar S-wave processes nn — >■ nn, KK, rjrj and on decays J/tp 
(pTTTT, 4>KK from the DM2, Mark3 and BES collaborations is performed for studying /o-mesons. 
The method of analysis is based on analyticity and unitarity and uses an uniformization procedure. 
In analysis only of the multi-channel TTTr-scattering, two possibilities were found for parameters of 
the /o(600) with mass about 700 MeV and with total width about 600 and 930 MeV. The use of 
only DM2 and Mark3 data on the J/ip decays does not permit to choose any of these possibilities. 
However the di-pion mass distribution in the decay J/ijj ^ 4>Tr'^n~ from BES prefers clearly the 
wider /o(600). Some spectroscopic implications from results of the analysis are discussed. 
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I. INTRODUCTION 

The comprehension of nature of /o mesons is very important for such profound topics in particle physics as the 
QCD vacuum. However both parameters of the scalar mesons, obtained from experimental data in various analyses, 
and even the status of some of the mesons, are not quite definite P, When discussing the first point, we pick out 
the /o(600)/(T- meson, /o(980) and /o(1500). As to the second point, one might indicate a situation related to the 
/o(1370) when, e.g., D. Bugg has indicated a number of data requiring apparently the existence of the /o(1370). 
These are the Crystal Barrel data on pp — ^ 777777° and on pp — >■ Stt" also the BES data on J/ip — ^ (/)7r+7r~, and in 
the GAMS data for tt+tt" — ?> tt^tt^ at large \t\. On the other hand, in works j^] one did not find evidence for the 
existence of the /o(1370). We have shown 2] that the existence of the /o(1370) does not contradict the data on 

TTTT — ^ TTTT KK, rjl], rjTj' . 

Note a situation with scalar states in the 1500- MeV region. In our previous mo del- independent analyses of tttt — > 
TTTT if if , 7777(7777') we saw the wide state /o(1500) whereas in works of some other authors, analyzing mainly mesons 
production and decay processes and cited in the PDG tables [ij, the rather narrow /o(1500) was obtained. We 
suggested that the wide /o(1500), observed in the multi-channel tttt scattering, indeed, is a superposition of two 
states, wide and narrow. The latter is observed just in decays and productions of mesons. Here we verify also this 
suggestion. 
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In view of indicated circumstances, related to parameters and status of scalar mesons, there are known problems 
as to determining their QCD nature and assignment to the quark-model configurations in spite of a big amount of 
work devoted these problems (see, e.g., and references therein). 

We present results of the combined coupled-channel analysis of data on isoscalar S-wave processes tttt — ^ tttt, KK, rjri 
and on decays J ftp — ^ tttt, KK for studying /g-mesons lying below 1.9 GeV. To analysis of experimental data we applied 
our model-independent method based only on the first principles (analyticity and unitarity) 6] . That approach permits 
us to omit theoretical prejudice in extracting the resonance parameters. Considering the obtained arrangement of 
resonance poles on the Riemann-surface sheets, obtained coupling constants with channels and resonance masses, we 
draw definite conclusions about nature of the investigated states. 

II. THE 3-COUPLED-CHANNEL FORMALISM IN MODEL-INDEPENDENT APPROACH WITH 

UNIFORMIZING VARIABLE 

Our model-independent method which essentially utilizes an uniformizing variable can be used only for the 2- 
channel case and under some conditions for the 3-channel one. Only in these cases we obtain a simple symmetric 
(easily interpreted) picture of the resonance poles and zeros of the S'-matrix on the uniformization plane. The 3- 
channel S'-matrix is determined on the 8-sheeted Riemann surface. The matrix elements Sij, where i,j = 1,2,3 
denote channels, have the right-hand cuts along the real axis of the s complex plane (s is the invariant total energy 
squared), starting with the channel thresholds Si [i — 1, 2, 3), and the left-hand cuts. The Riemann-surface sheets are 
numbered according to the signs of analytic continuations of the square roots — Si (i = 1, 2, 3) as follows: 

signs ^Im-^s — si, lm^/s^^J^, Im^/s — = + + +, — h+, l-,H , , — I — — correspond to 

sheets I, II, - • • , VIII, respectively. 

The resonance representations on the Riemann surface are obtained with the help of formulas from 0, @ , expressing 
analytic continuations of the S'-matrix elements to all sheets in terms of those on sheet I that have only the resonances 
zeros (beyond the real axis), at least, around the physical region: 

TABLE I: Analytic continuations of the 3-channel S-matrix elements to unphysical sheets 



Process 


1 ^ 1 


1^2 


2^2 


1 ^ 3 


2^3 


3^3 


I 


5ii 


S'12 


S22 


513 


S23 


<S'33 


II 


1/Sii 


i<S'i2/5'ii 


D33/S11 


iSi3/Sii 


D23/S11 


£'22/511 


III 


S22/-D33 


—S12/D33 


S11/D33 


-iDi3/D33 


iD23/D33 


det S/D33 


IV 


D33/S22 


iSl2/S22 


I/S22 


— -D13/S22 


iS23/S22 


D11/S22 


V 


det S/Dii 


iDi2/Dii 


533/1^11 




~S23/Dll 


S22/D11 


VI 


Dii/detS 


-Di2/detS D22/detS 


Di3/detS 


-D23/detS D33/detS 


VII 


S33/D22 


iDl2/D22 


det 5'/_D22 


— S13/D22 


iD23/D22 


S11/D22 


VIII 


D22/S33 


D12/S33 


Dii/5'33 


iSl3/S33 


*523/S'33 


l/5'33 



In Table HI the superscript / is omitted to simplify the notation, det S is the determinant of the 3x3 S- matrix on 
sheet I, Dai3 is the minor of the element Sa^, that is, Dn = S22S33 — S23, D22 = S11S33 — Sf^, D33 = S'iiS22 — S'12, 
D12 = S12S33 — S13S23, D23 = SiiS'23 — Si2S'i3, etc. 

Then, starting from the resonance zeros on sheet I, one can obtain an arrangement of poles and zeros of resonance 
on the whole Riemann surface. 

In the 3-channel case, we obtain 7 types of resonances corresponding to 7 possible situations when there are resonance 
zeros on sheet I only in Sn - (a); S22 - (b); S33 - (c); Sn and S22 - (d); S22 and S33 - (e); Sn and S33 - (f); 
Sii, S22, and S'33 - (g). 

The resonance of every type is represented by the pair of complex-conjugate clusters (of poles and zeros on the Rie- 
mann surface). A necessary and sufficient condition for existence of the multi- channel resonance is its representation 
by one of the types of pole clusters. Whereas cases (a), (b) and (c) can be related to the resonance representation by 
Breit-Wigner forms, cases (d), (e), (f) and (g) practically are lost at the Breit - Wigner description. 

The cluster type is related to the nature of state. E.g., if we consider the tttt, KK and rjrj channels, then a resonance, 
coupled relatively more strongly to the tttt channel than to the KK and rjrj ones is described by the cluster of type 
(a). In the opposite case, it is represented by the cluster of type (e) (say, the state with the dominant ss component). 
The glueball must be represented by the cluster of type (g) as a necessary condition for the ideal case. 
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From formulas of the analytic continuations, we conclude that parameters of resonances (masses, total widths and 
coupling constants with channels) must be calculated using the pole positions on sheets II, IV and VIII because the 
analytic continuations only onto these sheets have the forms: cx l/S\i, cx l/S'22 and oc 1/533, respectively, i.e., the 
pole positions of resonances only on these sheets are at the same points of the complex-energy plane, as the resonance 
zeros on the physical sheet, and are not shifted due to the coupling of channels. 

One can distinguish 0,0], in a mo del- independent way, a bound state of colourless particles {e.g., KK molecule) 
and a qq bound state. Just as in the 1-channel case, the existence of the particle bound-state means the presence of a 
pole on the real axis under the threshold on the physical sheet, so in the 2-channel case, the existence of the bound- 
state in channel 2 {KK molecule) that, however, can decay into channel 1 (tttt decay), would imply the presence of 
the pair of complex conjugate poles on sheet II under the second-channel threshold without the corresponding shifted 
pair of poles on sheet III. In the 3-channel case, the bound state in channel 3 (7777) that, however, can decay into 
channels 1 (tttt decay) and 2 {KK decay), is represented by the pair of complex conjugate poles on sheet II and by 
the pair of shifted poles on sheet III under the -qr] threshold without the corresponding poles on sheets VI and VII. 

It is convenient to use the Le Couteur- Newton relations 8]. They express the iS-matrix elements of all coupled 
processes in terms of the Jost matrix determinant d{^ s — si, • • • , ^/s — s„) that is a real analytic function with the 
only square-root branch-points at s ~ Si — 0. The important branch points, corresponding to the thresholds of 
the coupled channels and to the crossing ones, are taken into account in the proper uniformizing variable. Here we 
used an uniformizing variable, in which we neglect the lowest TTTr-threshold branch-point and take into account the 
threshold branch-points related to two remaining channels and the left-hand branch-point at s = 

v/(s-g2)53 + V(g- 53)52 
^S(S3 - S2) 

here S2 = S3 — Am'^. It maps our model of the 8-sheeted Riemann surface onto the uniformization w-plane. 

On Figs, im resonances of the (a), (b), (c) and (g) types, met in this analysis, are represented in 6*11 by the poles (*) 
and zeros (o) symmetric to these poles with respect to the imaginary axis. The Roman numerals denote images of 
the corresponding sheets of the Riemann surface; the thick line represents the physical region; the points tttt, i, 1 and 
t' = {\/^ + y/s2)/V^3 ~ 52 correspond to the tttt, KK, rjrj thresholds and s = 00, respectively; the shaded intervals 
(— CX),— 6], [—6^^, h~^,[b, 00) are the images of the corresponding edges of the left-hand cut of the TTTr-scattering 
amplitude. 

On the w-plane, the Le Couteur-Newton relations are [gI]: 

_ d*{-w*) _ d{-w-') _ d{w-') 

^''-^{wT' d{w) ' '^^^-^(^' 

Q Q q2 _ d*{w*~'^) 2 _ d*(-w*"^) 

»J11'J22 ^ •J12 — J? — \ ' 'Jll'-'33 ^ '-'13 



d{w) ' d{w) ■ 

The main model-independent effect of multi-channel resonances is given by the pole clusters. Assuming that possible 
remaining small (mo del- dependent) contributions of resonances can be included via the background, the S'-matrix 
elements are taken as the products 

S = SsSres (3) 

where Sb describes the background, Sres the resonance contributions. The d-function for resonance part is 

M 

dres{w) = W--^ Y[{W + w;) (4) 



with M the number of resonance zeros. For the background part Sb, the c?-function has the following form 

s — t 
2m„ 



d_B = exp[-i(a + ^ — -{an + if^n))] (5) 



n=l 

where 

s — s s — s 

an = flril + ana -d{s - Sa) + a„t, -9{s - Sv), (6) 

S ' — s s — S- 
Pn = bnl + hna -d{s - S^) + hnv -0{s - S„) 

with Sa the acr threshold, Sy the combined threshold of many opening channels in the range of ~ 1.5 GeV (????', pp, ujuj). 
These threshold are determined in the analysis. 
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FIG. 1: Uniformization w-plane: Representation of resonances of types (a), (b), (c) and (g) in the 3-channel vrTr-scattering 
5- matrix element. 



III. ANALYSIS OF THE DATA ON ISOSCALAR S-WAVE PROCESSES tttt -> tttt, KK, rjrj AND ON 

DECAYS J/ip (ji-RTx, <t>KK 

In the combined analysis of data on processes tttt — > mr^KK, rjri (references to used data are in Q) we added also 
data on decays J/ifi (pinr, (pKK from Mark III DM2 and BE S fill collaborations. Formalism for calculating 
di-meson mass distributions of these decays can be found in Refs. 0, Il2l|. There it is assumed that the pairs of the 
pseudo-scalar mesons of the final states have / = J = and that only these pairs undergo strong interactions, whereas 
the (f> meson acts as a spectator. 

The amplitudes for J /ip — >■ ipTni, (pKK decays are related with the scattering amplitudes i,j — 1 — tttt, 2 — KK 
as follows 

F{J/^ ^ 0w) = [ci(s)Tn + C2(s)r2i], (7) 

FiJ/^^ 4>KK) = v/l72 [ci (s)Ti2 + C2 (s)r22] , (8) 

where q = OLi/[s — + 7^0 + lus are functions of couplings of the J/il) to channel i; a^, /3i, 7^0 and 7^1 are free 
parameters. The first term in Ci is related to the on-shell manifestation of Adler zero. The di-meson mass distributions 
are given as 

N\F\''^{s - .OK - - ^1>f)i^l - + ^t>?) (9) 

where N (normalization to experiment) is 0.725 for Mark III, 0.347 for DM2 and 5.653 for BES data (from analy- 
sis). Parameters of the Q-functions, obtained in the analysis, are ai,/3ia2,/32 = 0.0306, 0.0646, 0.0222, 0.0701 and 
710,711,720,721 = 1.1902, 1.3622,-1.2823,-1.6748. 

In the analysis we supposed that in the 1500-MeV region there are two resonances: the narrow /o(1500) and wide 
/g(1500). From the analysis, the /o(600) is described by the cluster of type (a); /o(1500), type (c) and /o(1500), type 
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(g); the /o(980) is represented only by the pole on sheet II and shifted pole on sheet III. However, for representation 
of the /o(1370) and /o(1710), there are some possibilities. These states are described by clusters either of type (b) or 
type (c). Analyzing only the processes tttt — >■ tttt, KK, rjrj, rjrj' it is impossible to prefer surely any of four indicated 
possibilities; moreover, it was found that the data admit two possibilities for parameters of the /o(600) with mass, 
relatively near to the p-meson mass, and with the total widths about 600 and 950 MeV. 

When adding to the combined analysis, the data on decays J/V' — (/)7r7r, (pKK, one can give some preference for 
scenarios when the /o(1370) is described by the cluster of type (b), /o(1710) either of type (b) or (c). Further, for 
definiteness, we shall tell about the case when the /o(1710) is represented by the cluster of type (c). 

It is interesting that the di-pion mass distribution of the J/tp (piriT decay of the BES data from the threshold 
to about 0.85 GeV prefers surely the solution with the more wide /o(600). Satisfactory combined description of 
all analyzed processes is obtained with the total x^/NDF = 424.317/(389 - 55) « 1.26; for the TTTT scattering, 
xVNDF w 1.20; for vrvr KK, xV^DF « 1.64; for tttt rj-q, xV^dp ~ 0.82; for decays J/'0 ^ (j)T:iT , (j)KK , 
xVNDF « 1.55. 

Figures 2-6 demonstrate the fitting to experimental data. 
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FIG. 2: The phase shift and module of the vrvr-scattering S-wave matrix element. 
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FIG. 3: The phase shift and module of the tttt — > KK S-wave matrix element. 
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FIG. 4: The squared modules of the tttt — >■ rj-q S-wave matrix element. 



In Table HIl we show the obtained pole clusters for the resonances in the complex energy plane y^s; the poles on 
sheets IV, VI, VIII and V, corresponding to the /o(1500), are of the 2nd and 3rd order, respectively (this is an 
approximation) . 
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FIG. 5: The J/tp ^ ^tttt, (fiKK decays. The upper panel shows the fit to data of Mark III, the lower to DM2. 
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FIG. 6: The J/tp — > (^tttt decay; the data of BES Collaboration. 



The obtained background parameters are: a = an — 0.0, ai^ = 0.0198, aiy = 0.0, bn = bi^ ~ 0.0, biy — 0.0336, 
= -2.4808, a2a = -2.3021, a2v = -6.62, &21 = 62,7 = 0.0, b2v = 6.99, 631 = 0.6432, 63^ = 0.489, &2t, = 0; 
s„ = 1.6384 GeV^ s„ = 2.0851 GeV^. 

The obtained very simple description of the TTTr-scattering background confirms well our assumption (j3|). 

Generally, wide multi-channel states are most adequately represented by pole clusters, because the pole clusters 
give a main model-independent effect of resonances. The pole positions are rather stable characteristics for various 
models, whereas masses and widths are very model-dependent for wide resonances. However, mass values are needed 
in some case, e.g., in mass relations for multiplets. Therefore, we stress that they should be calculated using the poles 
on sheets II, IV, VIII, depending on the resonance classification. E.g., if 
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TABLE II: The pole clusters for resonances on the -^5- plane. y^=Er — irr/2. 



Sheet 




/o(600) 


/o(980) 


/o(1370) 


/o(1500) 


/o(1500) 


/o(1710) 


1— 1 
1— 1 


E, 

r./2 


521.6 ± 12.4 
467.3 ±5.9 


1008.4 ±3.1 
33.5 ±1.5 






1512.4 ±4.9 
287.2 ± 12.9 




III 


r./2 


552.5 ± 17.7 
467.3 ±5.9 


976.7 ±5.8 
53.2 ± 2.6 


1387.2 ± 24.4 
167.2 ±41.8 




1506.1 ±9.0 
127.8 ± 10.6 




IV 


E, 

r./2 






1387.2±24.4 
178.2 ±37.2 




1512.4±4.9 
215.0 ± 17.6 




V 


E, 

r./2 






1387.2±24.4 
261.0 ±73.7 


1493.9 ±3.1 
72.8 ±3.9 


1498.8 ± 7.2 
142.3 ±6.0 


1732.8 ±43.2 
114.8 ±61.5 


VI 


r./2 


573.4 ± 29.1 
467.3 ±5.9 




1387.2±24.4 
250.0 ±83.1 


1493.9 ± 5.6 
58.4 ± 2.8 


1511.5 ±4.3 
179.3 ±4.0 


1732.8±43.2 
111.2 ±8.8 


VII 


Er 

r./2 


542.5 ±25.5 
467.3 ±5.9 






1493.9 ± 5.0 
47.8 ± 9.3 


1500.4 ±9.3 
99.9 ± 18.0 


1732.8±43.2 
55.2 ±38.0 


VIII 


Er 

r,/2 








1493.9 ±3.2 
62.2 ±9.2 


1512.4±4.9 
298.4 ± 14.5 


1732.8±43.2 
58.8 ± 16.4 



TUres = \J + i^r/'^f' and Ttot = Tr- 

In Table Hill we show the obtained masses and widths for the /o states. 



TABLE III: Masses and total widths of the /o states. 





/o(600) 


/o(980) 


/o(1370) 


/o(1500) 


/o(1500) 


/o(1710) 


mres[MeV] 


700.3±10.0 


1009.0±3.1 


1398.6±24.7 


1495. 2±3. 2 


1539.4±5.4 


1733.8±43.2 


rtot[MeV] 


934.6±11.8 


67.0±3.0 


356.4±74.4 


124.4±18.4 


574.4±25.8 


117.6±32.8 



IV. DISCUSSION AND CONCLUSIONS 



• In the combined model-independent analysis of data on tttt — > tttt, KK^ rj-q in the I'~^J^'~" — 0+0++ channel and on 
J/ip ^ (pTTTT, (pKK from Mark III, DM2 and BES collaborations, an additional confirmation of the /o(600) with 
mass about 700 MeV and width 930 MeV is obtained. This mass value accords with prediction (too- « uip) on the 
basis of mended symmetry by Weinberg [31 ■ This is also in agreement with a refined analysis using the large- -/Vj, 
consistency conditions between the unitarization and resonance saturation suggesting — m„ ~ 0{N^^) pTs! ]. 

• Indication for /o(980) {nires = 1009 MeV, Ttot = 67 MeV) is obtained to be, e.g., the bound r\r\ state, because 
this state lies slightly above the KK threshold and is described by the pole on sheet II and by the shifted pole 
on sheet III without the corresponding (for standard clusters) poles on sheets VI and VII. 

• The /o(1370) and /o(1710) have the dominant ss component. Conclusion about the /o(1370) quite agrees with 
the one of work of Crystal Barrel Collaboration [lB| where the /o(1370) is identified as r\r\ resonance in the tt^t??? 
final state of the annihilation at rest. This explains also quite well why one did not find this state considering 
only the tttt scattering Conclusion about the /o(1710) is consistent with the experimental facts that this 
state is observed in 77 — > KsKs and not observed in 77 — > tt+tt^ [l^. 

• In the 1500-MeV region, indeed, there are two states: the /o(1500) (m^es ~ 1495 MeV, Ttot ~ 124 MeV) and 
the /o(1500) {rrires ~ 1539 MeV, Ttot ~ 574 MeV). The /o(1500) is interpreted as a glueball taking into account 
its biggest width among enclosing states [l^. 
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• We propose the following assignment of the scalar mesons to lower nonets, when excluding the /o(980) as the 
non-qq state. The lowest nonet: the isovector ao(980), the isodoublet Kq{800), and /o(600) and /o(1370) as 
mixtures of the 8th component of octet and the SU(3) singlet. The Gell-Mann-Okubo (GM-0) formula 

3m% = ^m\. - ml^ 

gives m/g = 870 MeV. 

In relation for masses of nonet 

the left-hand side is by about 17% bigger than the right-hand one. 

• For the next nonet we find: the isovector ao(1450), the isodoublet i^Q(1450), and two isoscalars /o(1500) and 
/o(1710). From the GM-0 formula, ruf^ « 1450 MeV. In formula 

™/o(1500) +TO/o(i710) = 2mifj(i45o) 

the left-hand side is by about 11% bigger than the right-hand one. 

• This assignment removes a number of questions, stood earlier, and does not put any new. The mass formulas 
indicate to non-simple mixing scheme. The breaking of 2nd mass relations tells us that the cr — /o(1370) and 
/o(1500) — /o(1710) systems get additional contributions absent in the Xq(900) and i4rQ(1450), respectively. A 
search of the adequate mixing scheme is complicated by the fact that here there is also a remainder chiral 
symmetry, though, on the other hand, this permits one to predict correctly, e.g.^ the cr-meson mass [14] . 
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